Isospin Breaking in K£4 Decays of the 
Neutral Kaon 

A. Nehme 
February 1, 2008 



Centre de Physique Theorique 
Luminy, case 907 
F- 13288 Marseille Cedex 09 

nehmeOcpt . univ-mrs . f r 

Abstract 

In the presence of photons, the neutral Kg^ decay, — > 7r'^7r~£+f^, 
can be parameterized in terms of three vectorial, one anomalous, and 
one tensorial form factors. We present here analytic expressions of 
two vectorial form factors, / and g, calculated at one-loop level in 
the framework of chiral perturbation theory based on the effective La- 
grangian including mesons, photons, and leptons. These expressions 
may then be used to disentangle the Isospin breaking part from the 
measured form factors and hence improve the accuracy in the deter- 
mination of vrvr scattering parameters from K£4^ experiments. 

keywords: Electromagnetic Corrections, Kaon Semileptonic Decay, Form 
Factors, Chiral Perturbation Theory. 
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1 Introduction 

The fundamental state for Quantum Chromodynamics (QCD) remains a 
poorly understood medium due to its complicated non perturbative struc- 
ture. It is widely believed that one of the main features of the vacuum 
structure is color confinement. In the chiral limit, defined by the vanishing 
of the light quark masses, confinement leads to the spontaneous breaking of 
chiral symmetry provided that the number of light flavors, Nf, is at least 
three pQ. Assuming further that the QCD vacuum angle vanishes, then, the 
fundamental state is invariant under the vector (residual) symmetry ,2]. As 
stated by Goldstone's theorem Nj — 1 massless pseudoscalars, 0°, are 
created from the vacuum and coupled to the axial currents, A'^, through 



The pseudoscalars form a non linear realization of chiral symmetry j3] which 
becomes linear when restricted to vector transformations 'S]. Furthermore, 
their dynamics is dictated by the present symmetry breaking pattern which 
constrains them to interact via derivative couplings [H]. Relaxing the hy- 
pothesis of vanishing quark masses in favor of the assumption that quarks 
are massive but light, the pseudoscalars acquire dynamical masses by means 
of the vacuum alignment mechanism [3 |H] . This constitutes the basis for a 
systematic low-energy expansion in powers of p/A and rUq/A where q stands 
for a light quark (m, d, and probably s), and A is the mass threshold for pro- 
ducing particles heavier than the pseudoscalars. Taking Nj = 3 and identify- 
ing the pseudoscalars with the eight light mesons of the hadronic spectrum, 
the preceding (chiral) expansion is nothing else than the effective field theory 
of QCD and is called Chiral Perturbation Theory (ChPT) PHU]. AH the 
features of the underlying theory, and particularly its symmetry breaking 
pattern, are coded in the coefficients of the chiral expansion which are the 
parameters of the effective Lagrangian. Thereby, they are called order param- 
eters of the chiral symmetry since they point out its spontaneous breaking. 
Among the infinitely many order parameters the coupling Fo and the quark 
condensate, (qq), are of particular interest. While the non vanishing of the 




(1) 
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former is necessary and sufficient for the spontaneous breaking of chiral sym- 
metry, a non zero value for the latter breaks chiral symmetry spontaneously 
but the converse is not necessarily true ^I]. Chiral symmetry constraints are 
unable to predict the values of order parameters. In return, large Nc count- 
ing rules plead in favor of a sizable quark condensate. In order to confirm or 
reject this argument, a phenomenological determination of order parameters 
is primordial. With respect to the coupling Fq, it can be extracted from the 
pion decay constant, 

F^ = Fo [1 + 0{m,)] . (2) 

Concerning the quark condensate, it enters the expression of all observables 
through meson masses [T^ . 

F^M^ = -(m„ + m,){qq) + 0{ml) . (3) 

Therefore, one cannot determine the size of this order parameter without a 
preliminary assumption on the chiral order of quark masses. The large Nc 
prediction means that the pion mass is dominated by the contribution from 
the quark condensate. Consequently, the first term in the right-hand-side of 
(ini) is leading order, that is, 0{mq) = O(p^). It follows that the ratio 

is proving a suitable parameter to test the present power counting scheme. 
Important deviations of X from 1 should incite on a careful reexamination 
of the validity of large arguments in the vacuum sector. It is interest- 
ing then to perform accurate measurements of observables that are sensitive 
to variations of X. This is the case of low-energy tttt scattering. In fact, 
the latter is solely described in the threshold region in terms of the 5'-wave 
scattering lengths, 

= [5« + 16/3 + 0K)] , (5) 



where, 



a 



4 - 3X + C(mq) , p = 1 + Oirrig) . (7) 
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For instance, the isoscalar scattering length, af], varies by a factor 2 when X 
varies from 1 to zero. This sensitivity of tttt scattering to the size of the quark 
condensate is not only a feature of the leading order but, on the contrary, it 
persists at higher orders allowing tttt scattering to be the golden process for 
testing the mechanism of quark condensation \VS\ . 

The TTTT scattering is experimentally accessible in K^^^ and pionium decays. 
In fact, the partial wave expansion of form factors displays the tttt phase 
shifts as stated by the Watson final state interaction theorem I15j. The 
data on phase shifts can then be translated into a model-independent deter- 
mination of scattering lengths by means of Roy equations PS]- Concerning 
the pionium characteristics such as its lifetime, r, and 2S — 2P strong energy 
level shift, AEg, they give direct access to vrvr scattering lengths via [T7] . 

r oc (a[! - alf , AEs oc 2a[! + al . (8) 

Once the results from the presently running DIRAC experiment are avail- 
able, Cq and Oq should be determined with 5 to 10% accuracy [TH]. On the 
other hand, the charged K^^ decay has been measured by the E865 experi- 
ment [ini 1201 outgoing data have been analyzed in |^ and [22] , in- 
dependently. Before recalling the conclusions of both references, let us stress 
that the obtained value for Og is 7% accuracy and is compatible with the 
prediction of the standard power counting which rests on large A^^^ grounds. 
The analysis of Ref. |2I] relies on chiral symmetry constraints which corre- 
late the two scattering lengths. It yields a value for Oq consistent with the 
standard counting prediction. Furthermore, if one combines the extracted 
values for the scattering lengths, the ratio Q reads then, X ~ 0.94. As 
for the analysis performed in Ref. it combines the data of the E865 
experiment with an existing one in the isospin-two channel below 800 MeV 
without using the aforementioned chiral constraints. The conclusions of this 
analysis point out a discrepancy at the 1 — a level between the extracted 
value for Oq and the one predicted by the standard power counting. More- 
over, the corresponding value for the ratio (jH) was found to be, X ~ 0.81. 
In view of the "disagreement" between the results of the two analysis, and 
before drawing off any conclusion about the size of the quark condensate, it 
is necessary to devote much more effort at both the experimental and the 
theoretical levels. In this direction, new precise measurements of charged and 
neutral Kqa, decay are currently taking data at CERN [22] and FNAL 
respectively. This should be accompanied by an improvement of the accuracy 
in the theoretical prediction for the scattering lengths by evaluating isospin 
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breaking effects in vrvr scattering as well as in decays. While such ef- 
fects in the former case are now under control at leading and next-to- leading 
orders [23 123 123 123 123 EQ|, we are interested in evaluating these effects 
at the same orders for the latter case. The aim of the present work is to 
evaluate isospin breaking effects in the Ki^ decay of the neutral kaon. The 
evaluation of the same effects in the decay of the charged kaon are now 
under control and will be published soon. 

2 decays in the isospin limit 

The semileptonic Kg4^ decays are given schematically by, 

K{p) 7r(pi) n{p2) t{pi) v,{pu) , (9) 

where the lepton, i, is either a muon, /x, or an electron, e, and u stands for 
the corresponding neutrino. 

Let us start with a historical survey on K14 decays calculation taking into 
account two important events (hence, three periods): The advent of Current 
Algebra and that of Effective Lagrangians. 

Before Current Algebra. The study of Ki^ decays starts in the late 
fifties jST] with an attempt to extract some information about the decay 
rate from comparison of the available volume in phase space and of the ma- 
trix elements for Ki4 and Ki^ decays. The interest was to test some models of 
the strong interaction. In this direction, the decay rate was computed on the 
basis of the universal Fermi interaction and some restrictions was imposed 
by strong interaction selection rules [321. The study of Ke^ decays proceeded 
with the calculation of the energy spectrum for leptons [23 , the energy spec- 
trum for pions [MJ , and the distribution over the angle between the directions 
of emission of the pion and the lepton [^3- The preceding studies were per- 
formed neglecting the interaction of pions in the final state. Deviations from 
this statement were considered in [^ and [HZ] using dispersion relation tech- 
niques. The results obtained in [23 showed that the shape of the momentum 
distribution is sensitive to the final state interaction of pions while the de- 
cay rate is affected slightly. On the other hand, it has been shown in [2Z] 
that, under reasonable assumptions, the Ku decay mass spectra for the two 
pions are completely determined by the partial S- and P-wave amplitudes 
for KK TiTx scattering. Hence, experimental data on Ki^, decays can be 
used to obtain some definite information on tttt and ttK interactions. To this 
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end, calculations were made again under the assumption that K^^ decays are 
mediated by pion-pion and pion-kaon resonances in [HHl EH] or sigma reso- 
nance in |1U]. Another approach was proposed in [^T] pointing out the 
usefulness of correlations between angular variables describing K14 decays 
in the determination of the difference between the S-wave and P-wave tttt 
scattering phase shifts. This approach was then detailed in 133] with a 
full kinematical description of K^^ decays. Finally, a complete review about 
Ke4 decay studies during this period can be found in jl^l liEj. 

The current Algebra period. The first Current Algebra calculation of Ki^ 
decay form factors was made in 07] for pion momenta put equal to zero. 
Different values were obtained for one form factor depending on which of 
the two pions in the final state was taken to be soft. The calculation was 
reconsidered by the author of jlH] US] who explained that the aforementioned 
difference is due to a missing piece in the previous calculation corresponding 
to a A' pole. The residue of the latter being determined from nK scattering 
and has been discussed in jJU] . The Current Algebra values for form factors 
have then been used to predict decay rates jSI]. In view of the Current Al- 
gebra prediction for the form factors, and with intent to test this prediction 
experiment aly, it is of considerable interest to determine form factors from 
experiment in a direct and economical way. It has been shown in |^ 13^ 
that both form factors and vrvr phase shift difference can be directly obtained 
from the measurement of the intensity spectrum under the assumption that 
5*- and P-wave dipion states dominate over the higher partial waves in Ki4^ 
decays. The decay rates calculated from the values of form factors predicted 
by Current Algebra showed some discrepancy with the experimentally mea- 
sured ones at that time. In order to account for this discrepancy in the Kqa 
case, the form factors were calculated again [53j by extrapolating decay am- 
plitudes from soft pion limits to the physical point using collinear dispersion 
relation techniques. The same method was applied later to the K^4 case 
in j^Sj. For a complete review about the use of Current Algebra to study 
KiA decays, we refer to 

The Effective Lagrangian period. The first application of the effective 
chiral Lagrangian method to the calculation of K^^^ form factors was done 
in [57 using the non-linear realization of chiral symmetry. The calculation 
of form factors with the help of an effective Lagrangian with linear realiza- 
tion of chiral symmetry was done in [SH] where the experimental values for 
decay rates has been reproduced thanks to the mixing of scalar qq mesons 
to gluonium. Within the framework of ChPT, the form factors were 
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calculated in jSHl EOl EI] to next-to-leading chiral order and in the isospin 
limit. A fit between the obtained results and the experimental values led to a 
determination of low-energy constants of the chiral Lagrangian under large- 
Nc assumption. A model-dependent calculation of form factors was carried 
out in jn21 using chiral Lagrangian with hidden local symmetry where vec- 
tor mesons, = 1~, play the role of gauge bosons. The results were in 
significantly worse agreement with experiment compared with ChPT results 
and suggested the inclusion of a scalr = 0"*" resonance in the nn channel. 
In the framework of generalized ChPT, where the ratio Q is a free param- 
eter, the Ke4 form factors were computed to leading order by the authors 
of [HSj who showed that high precision ^^"^4 experiments should allow for a 
direct measurement of the quark mass ratio 771^/(711^ + m^). In view of the 
increasing accuracy offered by the one-loop level evaluation of form factors, 
a reconsideration of the method proposed in [S2] to extract Tin phase shifts 
from K£4^ decay experiments is indispensable. This was achieved in [HI] with 
the conclusion that next-to-leading order corrections to the form factors in- 
duce a correction to the method not exceeding the 1%. The first estimation 
of two-loop corrections to the form factors was made by the authors 
of 1^ using dispersion relations. The same subject was also considered by 
the author of j^Sl with the help of the inverse- amplitude method. He used 
unitarity and dispersion relations together with the chiral expansion to inves- 
tigate the effects of two-loop corrections coming from tttt rescattering. The 
one-loop and partial two-loop results have then been used by the authors 
of [nn] who proposed another method to extract vrvr phases from Ki^ decay 
data. The method assumes a linear parametrization for the dependence of 
form factors on kinematical variables and allows a better exploitation of data 
than the one described in j52|,lb4j. Finally, the full two- loop level calculation 
of form factors has been done in [HZl IHHl ini] • 

2.1 Matrix elements 

The decays Q are described in terms of an invariant decay amplitude. A, 
defined via the matrix element, 

= I (27r)^<5(^) (pi +P2+Pi+Pu-V){-^A) , (10) 
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with the on-shell conditions, 



p2 = M\ pI = Ml, pI = Ml, pI = ml, pi = 0. (11) 

In the absence of electromagnetism, the decay © proceeds at leading order 
in Perturbation Theory through the exchange of a boson between two 
leptonic left-handed currents. For energies very small compared to My/±, 
such decays can then be described by the effective local Hamiltonian, 

G 

Tieff = —j= JI t^£7^(l ~ 7^)^ + hermitian conjugate . (12) 
v2 ^ 

In the preceding equation, 

J' = J]^.7ni-7')V^.,^,, (13) 

stands for the usual V — A weak current, Vij denotes the Cabibbo-Kobayashi- 
Maskawa flavor-mixing matrix elements, ipi represents a light quark flavor, u, 
d, or s, and finally, Gp is the so-called Fermi coupling constant. Let AQ, AS", 
and A/ denote respectively the change of charge, strangeness, and isospin of 
the current (fT!?jl . Since the latter satisfies the AS* = AQ rule, then, three 
physical modes for the decay © are possible jSSEBl, 

K+{p) r^^{pi)^~{p2)t{pi)u,{p,), (14) 

K^{p) ^ 7i\p,)Tx\p2)t{pi)vM, (15) 
K\p) Ti\pr)Ti-{p2)t{pi)u,{p,). (16) 

The corresponding decay amplitudes will be denoted by A^" , J^^, and , 
respectively. From now on, we will refer to the decay (jl4j) as the charged 
channel, the decay (fT3jl as the neutral channel, and the decay (fT^ as the 
mixed channel. 

2.2 Isospin decomposition 

Pions and kaons are isospin eigenstates with eigenvalues 1 and 1/2, respec- 
tively. It follows that the two pions in the final state of Q have a total 
isospin / = 0,1,2. Assuming that the current (fTSj) satisfies the A/ = 1/2 
rule, then, the two-pion system must be a state of isospin and 1 only. More- 
over, one can assimilate the current (jl3p to a spurion carrying an isospin 1/2 
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and think about the strong interaction part of K^^ decays as being a scat- 
tering between the spurion-kaon system and the pion-pion one. Accordingly, 
the decay amphtudes for the reactions (fT^ - (fT^ can be expressed in terms of 
the invariant amphtudes, and A^, for Ki^ transitions to states of definite 
isospin, 1 = and / = 1, respectively ^, 

A- = Ia^-^^A\ (17) 

= -^A\ (18) 

= -^AK (19) 

This means that the amplitudes for the various processes are related, 

A+- = (20) 

The Bose symmetry of the two-pion system reflects an invariance under the 
interchange of the two pions in the final state, pi P2- Therefore, the 
isospin wave function of the two-pion system is symmetric for 1 = and 
antisymmetric for 1 = 1 under the interchange, pi ^ p2- Moreover, in 
accordance with Bose principle, the orbital angular momentum of the two- 
pion system should be even for / = and odd for / = 1. 



2.3 Form factors 

To see the consequences of the foregoing on the dynamics of Ki^ decays, let 
us consider the charged process, (HH), and introduce the notations, 

P = p^+p^^ Q = p^-p^^ L = pi + p^, N = pi-p^. (21) 

The decay amplitude for the process in question can be written from the 
foregoing as follows, 

^(j9i)7r^(p2)|^7^(l - l')u\K+{p)) . (22) 

^We use Condon-Shortley phase conventions. 
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Due to the opposit relative intrinsic parities of the and tx^t[~ states, the 
matrix element of the vector current between these two states transforms 
as an axial vector, whereas that of the axial current transforms as a vector. 
Therefore, the hadronic matrix element appearing in the preceding equation 
possesses the following Lorentz decomposition, 

{'K^{pr)Ti-{p2)\sru\K+{p)) 

= -JJ^'^^'^^^-^P^'Q'^^^'^ (23) 
{n+{p^)n-{p,)\sYl'u\K+{p)) 

Mk± ^ ! ^ \ J 

The K^4^ form factors for the charged decay, H^^ , , and are 

analytic functions of three independent Lorentz invariants which we denote 
by, 

= P^, U = {p-piY, = {p-p2f ■ (25) 

They are made dimensionless by inserting the normalizations, and 
M~±, in (|23p and (j24p . respectively. The fact that we have used the charged 
kaon mass is a purely conventional matter and corresponds to the choice of 
defining the isospin limit in terms of charged masses. 

The interchange pi ^ p2 is equivalent to t,r ^ Ut^. It is then convenient 
to introduce the combinations, 

2/^ = I (Syr, tyr, "Utt) i / (Stt; Wjr, ^tt) , 

/ = H,F,G,R, (26) 

where the I^^{St^, It^^ m^) are defined in (j23|) and (|24j) . Similar notations hold 
for the neutral and mixed channels. From the foregoing, the decay amplitude 
A^^ is symmetric under ^ Mtt, whereas J\^~ is antisymmetric. Thus, the 
corresponding form factors satisfy, 

^00 ^ ^00 ^ ^00 ^ ^00 ^ ^0- ^ ^0- ^ ^0- ^ ^0- ^ (27) 

Since fl20|) is also satisfied by the form factors, one obtains the following 
relations between the non vanishing parts in form factors of processes fllSp 
and (Uni), 

{H_ ,F+,G^, = -iH_,F+,G_, i?+)+- , (28) 
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,F.,G+, = V2{H+,F^,G+, . (29) 

Hence, it is sufficient to study one of the processes (fT^ - (fT^ in the isospin 
hmit. 1 



3 Isospin breaking 

Isospin breaking affects the decay rate as follows. It generates corrections to 
the form factors F, G, R and H. It introduces a supplementary form factor 
T. It modifies the Dalitz plot and contributes to the intensity spectrum. 



3.1 Kinematics 

We follow the approach used in to study the kinematics of K^^ decays. 
The idea consists on looking at these decays as being two-body decays into 
a dipion of mass s^r and a dilepton of mass si = L^. The two systems 
subsequently decay in their own center-of-mass frames. To describe the decay 
distribution it is convenient to use, besides the invariant masses, s^^ and si, 
the angles 6*^, 9i, and as illustrated in Fig. [TJ 

All of the scalar products obtained from momenta, pi, p2, pi, and p^, or 
equivalently from, P, Q, L, and A^, can be expressed in terms of the five 
independent variables, Stt, si, O-jr, 6i, and 0. In the kaon rest frame, the Dirac 
components of momenta read, 

p° = M, (30) 
pi = p2 ^ p3 = , (31) 

= -^J(M' + - + - Ml) 

+\'/^iM\s^,si)\'/''{s^,MlMl)cose^] , (32) 

+ {M^ + - si)\'/\s^,MlMl)cos9^] , (33) 
pI = -^X'/'is^,MlMl)sme^, (34) 
pI = 0, (35) 
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-X'/\M\s^,si)X'/\s,,MlM^)cose^] , (36) 
-(M^ + s^- si)X'/^{s^, Ml Ml) cos^,] , (37) 

P2 = 0, (39) 

= 4^[(s« + ^?)(M' + Si-s.) 

+(s,-m^)A^/2(M2,s,,sOcos^,] , (40) 
^' ^ -^[(«^ + mnAV2(M^s.,sO 

+ (si - mf)(M2 - + si) cosOi] , (41) 

P° = ^ (5; - W) [M^ -s^ + si- \"\Ml s,, si) cosOi] , (44) 

pI = {si - mf) [AV2(M^ s,, si) - (M^ - + s,) cos^,] , (45) 

pI = i^i - ^1) sin cos , (46) 

pI = (s; - ml) sin 6*; sin . (47) 

In the preceding equations, our choice for the metric tensor is, 

[rj] = diag{l, -1,-1,-1}. 
Therefore, it is easy to obtain the scalar products, 

^ s^, (48) 
Q' = ^-s^Y\ (49) 
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P-Q = Ai2, 

P-L = l{M^-s,-si), 

P-N ^ ]^{M^ -s^-si)zi + {l-zi)Xcosei, 

Q-L = {M^ - s^- si)p^ + XY cos e^, 



Q.N = (1 - — cos Oi + (M^ -s^-si 



2 s- 

+^iXF cos + ^ (1 - ^Oi^ [ (M^ -s^-si 

cos COS — sin sin 9i cos ] , 

L-N = ml, 
e^^p^L'^N^^P'^Q'' = -{l-zi)XYy/^ism9^sm9ism(f>, 

with the following notations, 

= — , 

Sl 

The function, 

X{x, y, z) = + ip + z^ — 2xy — 2xz — 2yz , 
is the usual Kallen function. 

3.2 Phase space integral 

The starting point is the K14 differential phase space, 
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d^p^ d^P2 d^Pi d^Pj^ 



{27r)^2Ei {2t:Y2E2 {27r)^2Ei {27rf2\p^\ 

with particle energies, 



E,^y/pf + pl (64) 

This can be written as 



= J d^P j d^L{27r)^S^^\P + L-p) 



X 



Performing the integrations over the 5 functions separately, we obtain for the 
phase space, 

Jd^ ^ ^(27r)5y" ds^J dsi\'/\M\s^,si)<^^^^i,^, 
with the dipion phase space, 

= j^-^>."'Hs„MlMl)jd(cos0,), (65) 
and the dilepton phase space. 

Prom the foregoing, the phase space can be obtained by integrating the 
differential phase space, 

d^^ = M^N{s^,si)dsT,dsid{cos6^)d{cos9i)d(l) , 

with, 
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over the following range of variables, 






< 





< 




2n 





< 


ei 


< 




7T 





< 


6^ 


< 






2 


< 




< 


(M 




2{Ml + M|) 


< 




< 


(M 


- mf ■ 



(67) 



3.3 The decay rate 

Electromagnetism breaks the V — A structure of Ki^ decay amplitudes. By 
Lorentz covariance, these can be written in general as follows, 



GfV* 



where. 



Mk± 

h T 



The quantities /, g, r, and h, will be called the corrected Ki^ form factors 
since their isospin limits are nothing else than the Ki^ form factors, F, G, 
R, and H, respectively. The normalization factors are written in powers of 
the charged kaon mass by convention (this corresponds to our definition of 
the isospin limit as stated in Sec. El). Squaring the amplitude and summing 
over spins, we get, 

J2\A\' = 2Gl\Vus\'^Ms^,si,e^,ei,(f)), 

spins 

with the following expression for the intensity spectrum, 

+ \g\' [(Q . Lf - {Q ■ - Q\si - mf)] 

+ \r\ nil ~~ "^i ) 

- {{e.^.^L^N^PPQ'^f + {si - mf) [Q^X' + s^{Q ■ Lf 

+ si{P-Qf-2{P-Q){Q-L){P-L)]] 
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+ i9*h + gh*)^[{P- L) {Q .L){Q.N)-{P. N) {Q ■ Lf 



+ {rg + fg*) [{P ■L){Q.L)-{P- N){Q . N) - {si - m^)iP ■ Q)] 
+ {fr + /r*)mf [(P • L) - (P • N)] 

+ {rh + fh*) ^ [iQ ■ iV)(P ■ Lf -{Q-L){P. L){P . N) 

- sMQ ■ N) + mfs^iQ ■ L) - m]{P ■ Q){P ■ L) + si{P ■ Q){P ■ N)] 
+ {g*r + gr*)m][{Q-L)-{Q-N)] 

1 

+ siQ''{P ■ N) - mfQ''{P ■ L) + (P ■ Q){Q ■ L) - si{P ■ Q){Q ■ N)] 

- e,.,^L>^N^P^Q" ^ [M\rg - fg*) 

- (P ■ N){rh - fh*) - {Q ■ N){g*h - gh*) - mUr*h - rh*)] 
+ {fT* + Tn^[s,{Q.N) 

- s^{Q-L) + {P.Q){P-L)-{P.Q){P.N)] 
+ {gT* + Tg*)^[Q\P.L) 

- Q'(P-N)-{P-Q){Q-L) + {P-Q){Q-N)] 

+ irT* + Tt*) ^ [(P ■L){Q.N)-{P. N){Q ■ L)] 

+ {rT* -Tr*)'^e^,p,L'^N''PPQ'' 

- {hT* + Th*) ^ [2(P ■ g)(P ■ L){Q ■L)-iP-Q){P. L){Q ■ N) 

- Q\P ■ Lf + Q\P ■ L){P ■ N) - s^{Q ■ Lf 

+ s^{Q ■L){Q-N)-{P- Q){P ■N){Q-L)- si{P ■ Qf 
+ {P-Qf{L-N) + s^siQ^-s^Q\L-N)\ 

+ \T? [s.siQ^ - si{P ■ Qf - 2Q\P ■ Lf 

+ 4(P • Q){P ■L){Q-L)- 2s^{Q ■ Lf - s^Q^N^ + (P • QfN^ 

+ 2Q2(p . Nf + 2s^{Q ■ Nf - 4(P • Q)(P • N){Q ■ N)] . (68) 

The Ki4^ differential decay rate is defined by, 

dr = -l^d^y\A\^. 

2M ^ ' ' 

spins 
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It reads, in terms of the five independent variables, 
d^T = Gp\Vus\'^N{s^, si)j5{s^, si,9^,9i,(f))ds-^dsid{cos9^)d{cos6i)d(f) . (69) 

3.4 The effective Lagrangian 

In order to treat completely electromagnetic effects in Ki^ decays, not only 
the pseudoscalars but also the photon and the light leptons have to be 
included as dynamical degrees of freedom in an appropriate effective La- 
grangian jTOj. The starting point is QCD in the limit m„ = = = 0. 
The resulting chiral symmetry, G = SU{3)l x SU{3)ji, is spontaneously bro- 
ken to SU{3)v- The pseudoscalar mesons {n,K,r]) are nothing else than 
the corresponding Goldstone fields 0j = 1, . . . , 8) acting as coordinates of 
the coset space SU{3)l x SU{3)r/ SU{3)v- The transformation rules for the 
coset variables ul,r{4>) are 

UR{<j)) giiUjih{g,(j))~^ , 
9 = {9L,gR) e SUi3)LxSUi3)R, (70) 

where h{g, 0) is the nonlinear realization of G E]- 

As stated before, the photon field and the leptons i, ui {i = e, /i) have 
to be dynamical. Thus, they most be introduced in the covariant derivative, 

Uf, = i[ u'^jiidf, - ir^)uR - u{{d^ - il^)uL ] , (71) 

by adding appropriate terms to the usual external vector and axial-vector 
sources V^, Afj_. At the quark level, this procedure corresponds to the 
usual minimal coupling prescription in the case of electromagnetism, and 
to Cabibbo universality in the case of the charged weak currents, 

= v^-a,-eQrA, + J2ih,^iLQl + m:i,iQV), 

= + - eQlf . (72) 

The 3x3 matrices QJ^ are spur ion fields corresponding to electromag- 

netic and weak coupling, respectively. They transform as. 
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under the chiral group. In practical calculations, one identifies with the 
quark charge matrix 

/ 2/3 
g^"" = -1/3 I , (74) 
\ -1/3 

whereas the weak spurion is replaced by, 



Ql = -2^ Gi 



us 



I , (75) 
\ 



where is the Fermi coupling constant and Vud, 14s a-re Kobayashi-Maskawa 
matrix elements. 

In order to take into account electromagnetic mass difference between 
charged and neutral mesons, it is convenient to define the following electro- 
magnetic and weak sources, 

L = KQl ^L, Qr = u'rQr Ur (76) 

transforming as 

QT'^ ^ h{gA)QT^h{gA)-\ 

QT ^ h{g,<p)QTh{9Ar'. (77) 

With these building blocks, the lowest order effective Lagrangian takes 
the form 

- ^F^^F'''' + Y^[i{i^ + e4-mi)i + VeLt^iyeL], (78) 

where ( ) denotes the trace in three-dimensional ffavour space. 

The low-energy constant Fq appearing in the preceding formula is an order 
parameter for chiral symmetry since it testifies to its spontaneous breaking. 

It represents the pion decay constant in the chiral limit, m„ = m^^ = = 0, 
and in the absence of electroweak interactions. Explicit chiral symmetry 
breaking due to quark masses is included in, 

X+ = UrXUl + u[x^ur . 
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In practice, one makes the following substitution, 

X ^ 25oA<quark = 2^0 Uld I , (79) 




where Bq is an order parameter for chiral symmetry. It is related to the 
quark condensate in the chiral limit by, 

{qq) = -F^Bo. (80) 

The low-energy constant Zq expresses explicit chiral symmetry breaking by 
electromagnetism. It is given by the electromagnetic mass of the pion as we 
will see below. 

In the absence of electroweak interactions, that is, for = , a = 
0, ChPT is Isospin- invariant. In order to study Isospin breaking effects in 
ChPT processes, the usual chiral expansion in powers of p and ruq is no more 
sufficient. One must also expand matrix elements in powers of the isospin 
breaking parameters, — m„ and a. On the other hand, the best accuracy 
ever reached in strong interaction observable measurements does not exceed 
the 5% level. Thus, we consider that an expansion to orders C(mrf — m„) and 
0{a) is highly adequate for our purposes. Moreover, chiral expansion and 
Isospin breaking expansion have to be related in a consistent way in order 
to obtain reliable results. We adopt an expansion scheme where the Isospin 
breaking parameters are considered as quantities of order in the chiral 
counting, 

0(m, - m„) = Oia) = 0(m,) = 0(/) . (81) 

Therefore, tree level calculation corresponding to leading chiral order is car- 
acterized by chiral orders cited in (|81|) . Concerning one-loop level calculation 
which corresponds to next-to-leading chiral order, it is caracterized by the 
following chiral orders, 

0{p') , 0{ml) , 0{p'm,) , 

O {p^md - m„)) , O (m,(mrf - m„)) , 0{p^a) , 0{mga) . (82) 

We have now all necessary elements to calculate any Green function in 
the framework of ChPT including Isospin breaking effects. For example, the 
leading chiral order expressions for light meson masses are found to be, 

M^o = Ml, (83) 
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(84) 
(85) 



(86) 




i (4M^ - Ml) 



(87) 



Herein, M^r and Mk represent respectively pion and kaon masses in the 
absence of isospin breaking, 

Ml = 2Bom, M|. = 5o(m + m,) , 2m = m„ + , (88) 

e measures the rate of SU (2) to SU (3) breaking. 



At next-to-leading chiral order, one-loop calculation is involved. Vertices are 
extracted from Lagrangian (fTSj) . Meson masses in the propagators as well 
as in the vertices can be identified with expressions (|HH j) - (jH7|) . As is well 
known, loops are ultraviolet divergent. To remove divergences, renormaliza- 
tion should be employed. The procedure consists on adding to Lagrangian 
(I7S|) suitable counter-terms [7TI [721 ESI EH EO] generating exactly the same 
divergences but with opposite sign. Moreover, the cancellation should occur 
order by order in the chiral expansion as dictated by renormalizability prin- 
ciples of effective field theories. Counter-terms are modulated by low-energy 
constants which are order parameters for chiral symmetry. In order to deter- 
mine these constants one proceeds as follows. Let C be either a low-energy 
constant or a combination of low-energy constants and T an observable very 
sensitive to variations of C. One first calculate the expression of T in the 
framework of ChPT to any given order and then match the obtained ex- 
pression with an experimental measurement of T. It is clear that the value 
for C deduced from this matching does not constitute a genuine determina- 
tion of the low-energy constant. In fact, it represents the value of C at the 
given chiral order and with the accuracy of the experimental measurement. 
Note that the method of effective Lagrangian has the disadvantage of an in- 
finitely increasing number of low-energy constants when going to higher and 
higher orders in the low-energy expansion. For instance, two constants in the 
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strong sector, Bq and Fq, and one constant in the electroweak meson sector, 
Zq, parametrize the leading chiral order. At next-to-leading chiral order, one 
has ten low-energy constants in the strong sector, Li , . . . , Lio, fourteen con- 
stants in the electroweak meson sector, Ki , . . . , K^, and seven constants in 
the electroweak leptonic sector, Xi , . . . ,Xi. The constants, Li, Ki and Xi, 
are divergent. They absorb the divergence of loops via the renormalization, 

Li = L[(/i) + r,A, ^ = 1,... ,10, (90) 
Ki = K[(/x) + S,A, t = 1,... ,14, (91) 
X, = X[(/i) + H,A, z = 1,... ,7. (92) 

Herein, A corresponds to pole subtraction in the MS dimensional regulariza- 
tion scheme (see appendix). The beta-functions, Fj, Sj and Hj, can be found 
in [7T], [72] and [7^, respectively. The scale cancels in observables as can 
be seen from the renormalization group equations, 

Llif^2) = LK/ii) + T^ln^, (93) 

and similar for Ki and X^. 



4 Leading order 

The charged, neutral, and mixed channels will be studied at tree level in- 
cluding isospin breaking effects. Due to the absence of a Lagrangian treating 
isospin breaking in the anomaly sector ^, the corrected form factor h will be 
put to its isospin limit, H. In the following, we shall give a detailed calcu- 
lation for the charged channel due to the presence of a tensorial form factor 
already at leading order. For the remaining channels, we limit ourselves to 
quoting results. 

4.1 Form factors 

The different topologies for K^ decays at tree level in perturbation theory 
are drawn in Fig. El 

The contribution of each Feynman diagram to the charged decay ampli- 
tude will be given separately. 

^See [21] for the two-flavor case. 
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diagram 2. (a) 

= -^G^C^(Pj(l + 7')7M(3i^'^ + 3g'^ + 2LXp,). 
diagram 2. (b) 

[ 3(s, + U- u^) - {si -Ml^) + 2AZoe^F^ ] L^v{pi) . (94) 
diagram 2. (c) 

diagram 2. (d) 

diagram 2. (e) 

= ^TO(P.)(l + 7-)^ ^.^_4_,^.^^ x 

+ (Q • L + Q • A^)7,.i^'' - 2imi(7^,p^p^] . 

Note the cancellation between the contributions from diagrams 2. 
(c) and 2. (e). Adding all this together, the corrected form factors in the 
charged channel read, 

r- = F^~, 



y/2Fo Si - 

2e^F^ f U - ,1^ Q- L + Q- N \ Mk± 



si - M^± M|,± -mf-2p-pj y/2Fo ' 



T = 



4e^Fo2 miMK± Mk± 



Sn M]^± -mf-2p-p^ y/2Fo ' 
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with the charged form factors, 



R+- 



Mk± 

Mk± 
2V2F0 
0. 



1 + 



si-Ml 



The corrected form factors in the neutral channel are given by, 

6e Mk± 

e Mk± 

1 











^00 





X 



1 + 



Si - M|± 



(3s^ + 2si - 6M^) 



with the neutral form factors, 

,00 



V2F0 



6-°° = 0, 

^00 ^ ^/^± 



2V2Fo 



if 



00 



0. 



The ones relative to the mixed channel have the following form, 

3e Mk± 

Mk± 







9'- 









e Mk± 



1 + 



1 



24 



(3s^ + 2si - 6M|-) 



with the mixed form factors 



F^' = 



G 



0- 



Mk± 
Fo ' 

o_ _ Mk± U - 
= 0. 



5 Mixed channel at next-to-leading order 

We present here a one-loop calculation of the Ki4 decay amplitude for the 
mixed channel including isospin breaking terms. Feynman diagrams repre- 
senting the amplitude will be separated into two sets: photonic and non 
photonic diagrams. The non photonic set is drawn in figure El 

The calculation of these diagrams is standard in field theory. The starting 
point is Lagrangian (ffF?|) . For the non linear realization of chiral symmetry 
we will use the exponential parametrization, 

ur = u[ = exp , (95) 



where $ is the linear realization of SU{3) and can be decomposed in the 
basis of Gellmann-Low matrices as, 



a=l 



In terms of physical fields the matrix $ can be written, 



$12 = -V2TT+ , $13 = , $21 = V2n 

$22 = (-l + ^]^'+(il + ^=]v, 



$23 = -V2K\ $31 = V2K-, $32 = -V^°, 

$33 = --^(e27r° + r/) . (96) 
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The two mixing angles ei and 62 relate 0^ and 0^ to the mass eigenstates 7r° 
and r], 

rO \ / 1 ^ \ / A3 



TT^ \ / 1 ei 

^ y V -e~2 1 



(97) 



Notice that ei = £2 = e at leading order. The next step consists on expand- 
ing Lagrangian (f75j) to fifth order in pseudoscalar fields, generating Feyn- 
man rules and drawing allowed topologies. We then calculate pseudoscalar 
propagators and derive masses and wave function renormalization constants. 
Finally, we expand the next-to-leading order Lagrangian to third order in 
pseudoscalar fields and obtain the counterterm contribution. 

Let us denote by SF and 6G the next-to-leading order corrections to the 
F^^ and G^~^ form factors, respectively, 

0- Mk± 



The expressions for 5F and 5G are lengthy. Therefore we will separate them 
to different contributions depending on the topology of the Feynman diagram 
representing a given contribution. 



5.1 Born contribution 

This contribution is obtained from diagram (a) in figure El We take the 
corresponding vertex from Lagrangian (f7S)) and multiply by the wave function 
renormalization constant factor to obtain, 



5F 



SG 



3e 

7! 



[19Ao(M^) + 3Ao(M,) + UAo{Mk)] 



F2 
-^0 



[3(M,2 + 2M|)L4 + {Ml + 2M'JL,] 



24F2 



6e \ 



Uo(M,o)+3 1 + 



2e 
73 



+ UAo{M^±) + 8Ao{Mko) + QAo{Mk±)] 



^ [3{M'^ + 2Ml)L, + (Mio + 2Ml,)U] 
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+ 2 



2 2 N 

5 + 21n— ^ + 21n^ 



- - {24Ki + 24^2 - 12X3 + 6X4 + I6X5 + 16^6 + SXe) . 

The Born contribution is infrared divergent. This divergence emerges from 
the wave function renormahzation constants of charged pion and lepton. 

5.2 ir^ - T] mixing contribution 

Isospin breaking induces ~ rj mixing by assigning non vanishing values to 
the off-diagonal matrix elements, 



(98) 



The calculation of these matrix elements is staightforward and one obtains 
for the two mixing angles the following next-to-leading order expressions, 



ei = e2 + 2e(Mi-M2)Ci 



3v^ (47r)2 



9Zo(l + ln^ 



62 



= e [1 - 3Ai^ 2/xk + Hr, 



32 



+ M^Cr ~ {Mf, - M^){3Lr + LI) 



Ml 



3v^ (47r)2 M^-Ml 



3Zo 1 + In 



Ml 



+ (471)2(67^3^ - 3KI - 2KI - 2KI + 2KI + 2K{^)\ 
Herein, the quantity Ci is defined by. 



1 



167r2F2 



1 



Ml 



Ml -Ml 

K TV 



1„M 

Ml 

7r 



and the tadpole integrals are denoted by, 

1 

lip 



(99) 
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Let us relate the two sets of mixing angles, (61,62) and (61,62). To this end, 
one calculates the pseudoscalar mass matrix to next-to-leading order using 
the parametrization and (jHSJ- Then the diagonalization condition leads 
to. 



61 



+ 



61 - ^ (Mi - Ml)C, 



1 + ln^ 



62 



(47r)2 V /U- 
e2 + ^ (Mi - Ml)C_ 



1 + ln^ 



(47r)2 



73 



Since 6 represents the — rj mixing angle at leading chiral order it is inter- 
esting to "renormalize" 6 by taking the mean value of the mixing angles and 
expand it to next-to-leading chiral order. 



2 i^i + ^2) 



+ 6 



(4) 



C(/ 



(100) 



where. 



6(^) 

^(4) 
^str. 



.(4) 



= 6 + + /i,, 

+ MlC,-—,{Ml-Ml){-,L, + U^ 
^0 



+ 



Mi 



9^3 (47r)2 M2-M2 [ 
2(47r)2(6ir3^ 
2e2 Ml 



9Zo 1 + In ■ 



M 



9^3 M2 - Ml 

[QKi - mi 



2KI 



2K: + QKl + QK 



10) 



(101) 



Note that the preceding equations were derived in [TH] neglecting terms pro- 
portional to e^M^. From the foregoing, one can treat 7r° — 77 mixing following 
two methods. The first consists on working with mass eigenstates, that is. 
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with pseudoscalars diagonalizing the mass matrix [77j. In practice, this is 
reached by keeping ii and 62 when deriving Feynman rules from Lagrangian 
fl7S|) with parametrization and (j^ . At the end one replaces ei and €2 
with their expressions quoted before. Notice that ei and €2 are divergent 
quantities and, thanks to this divergence, the final result is finite. We will 
not follow the method just described and opt for the following one. Let us set 
ii and €2 in parametrization (j95|) and (j96j) to their leading order value, e. As 
usual, derive Feynman rules from Lagrangian ()78|) and calculate one-particle- 
irreducible diagrams. After taking into account all contributions especially 
those coming from counter-terms, the final result still divergent naturally. 
This divergence can be cancelled by the one generated from diagram (b) in 
figure 121 which accounts for tt^ — rj mixing. The contribution of this diagram 
to the F and G form factors is, 



SF = 



+ 



{71 ^^^^ " ^^^'^ '^"^^'^"^ '71^^''' "^'^^'^ 

2M^o + M^o - ^ (Mi + 2M'^) 

2Mio + M,2„ - (Mi - 2M2) 
384e 



Ao(Mxo) 
Ao{Mk±: 



[Ml - Ml)\?,L, + L, 



- 2e^FlMl (6K3 - 3/^4 - 2K^ - 2K^ + 2Kc, + 2Ki^) ] 
5G = 



(102) 



We have checked that the two methods are completly equivalent. 



5.3 Counter-terms contribution 



The contribution in question follows from diagram (a) in figure El with the 
help of the next-to-leading order Lagrangian and reads. 



P-Pl-P-P2 + -^{P-Pl+P-P2-4:PI-P2) 



-2 4 [4{M' + '^MDU + 2(Mi + 2Ml)L, + s^L,] 
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+ - {-QKs + 3K4 + 2K5 + 2Ke - 6X1) , 



SG 



^0 L 



P ■ Pi + P ■ P2 + —7={P ■ Pi - P ■ P2) 



+ -F2 [4(M2 + 2M^)L4 + 2(M|o + 2M^2o)L5 + siL^] 
^0 

+ — {2AKi + 24ir2 - I8K3 + 9K4 + 2OK5 + 20^6 + ^^Ki2 - 6X1) . 
9 

5.4 Tadpole contribution 

Tadpoles are shawn in diagram (c) of figure El and contribute to the form 
factors by the following, 



5F 



1 



5e 



6e 



1 + -= Ao(M,o) + — Ao(M, 



V3 

Ao{M^±) + -^Ao{Mk±) 



V3 



SG 



1 



AFi 



2 1 



2e \ 



Ao(iV4o) + 1 - 



3e 

7! 



Ao(M,) 



+ 2 1- 



73 



Ao(M^o) + ( 3 + ) Ao(M,±) + 2Ao(M^±) 



5.5 The s-channel contribution 

The remaining contribution to the non photonic part of the decay ampli- 
tude comes from loop diagrams with two pseudoscalar propagators. This 
two-point function contribution will be separated in three parts depending 
on the Lorentz scalar governing its underlying kinematics. The s-channel 
contribution comes from diagram (d) in figure El and reads, 



6F 



3F2 



+ 



-AoiM^±)--^Ao{MK) 
3e 



M^o -2p,-p2 + -j= {Ml - 2pi ■ P2) 



Bo{-pi - p2,M^o,M^±] 



6e 

7! 



(M^ - Ml - pi ■ P2) Bo{-pi - P2, Mr,, M^) 
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+ 2 [SM^o - - Pi ■ p2 
3e 



6e 



5i(-pi-P2,M,o,M,±; 



6e 



+ -^{M],-MI-pi-P2)Bi{-pi~P2,Mk,Mk) 

+ 25oo(-pi - P2, M^o, M,±) + 5oo(-pi - P2, Mk, Mk) 
+ 4(2M2o - Ml± + pi ■ P2)5n(-pi - p2, M,o, M,±) 



- 3 



2e 



-P2,Mk,Mk) 



5G = -^[2Eoo(-pi-p2,M,o,M,±) + 5oo(-pi-p2,Mj^o,M^±)] 



5.6 The ^-channel contribution 

Diagram (f) in figure El generates the somewhat lengthy t-channel contribu- 
tion, 



SF 



1 



12Fi 



1 - ^] Ao{Mko) -3(1 + -=) Ao{Mk± 



V3j 



+ [2Mlo+p-pi 



3e 



+ —{4Ml-2Ml-p.p,) 

+ l[Ml-Ml,-2p-p^ 

+ ^{M^,-Ml + 2p.p,) 

- [Mlo + 2p-p^ 

+ ^(3M^-4M2-4p-pi) 

+ ^ [M^ + 2Ml, - SM^o 

+ ^[Ml-Ml + ml-2p.pr) 



Vs. 

Bo{pi-p,M^o,Mko) 
Boipi - p, M^, Mko) 



5i(pi-p,M,o,M^o) 



Bi{pi-p,M^,Mko) 
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+ 6 



+ 



5G 



+ 



+ 



p-pi--j^{2Ml-2Ml-p-pi) 



Bi{pi-p,M^±,Mk±) 



9e 



+ 4(1--^) BooiPl - p, M^o, Mko) 



+ 12 ( 1 - — J BooiPl - p, M^, Mko) 
+ 12 ( 1 - BooiPl - p, M^±, Mk±) 



+ [M^o - 2Mlo + p ■ pi 

- _(Mi-2M^+p-pi) 
+ 3 [M^o - 2M^o + p ■ Pi 

- ^{Ml-2Ml+p-pi) 

- &[p-pi-Mlo 



V3 



iMl-2Ml+p-pi) 



Biiipi-p,M^o,MKo) 



Biiipi-p,Mr„MKo) 



Biiipi-p,M^±,MK±) } , 



12Fi 



3e 



V3 

{AMl-2Ml-p-pi) 



Ao(M^o) + 3(l + — )Ao(M;^±) 



[2Ml,+p-pi 
3e 
73 

^-[Ml-Ml,-2p-pi 
3e 



Boipi-p,M^o,MKo) 



+ -^{M^,-Ml + 2p.pi) 

+ [M|.o + 2p • pi 
3e 

+ 



(3Mi - AMI - 4p • Pi) 



Biipi-p,M^o,MKo) 



- - [M^2 + 2Mio - 5M^o 
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{Ml-Ml + ml-2p-pi) 



Biipi-p,M„,MKo) 
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- 6 
+ 2|^1 
+ 6(1 



p-pi--^(2Mi-2M2-p.pi) 



9e 
e 

[M^o - 2M^o +P-PI 
9e 



5oo(Pi - P,M^o,Mko) 



Bi{pi-p,M^±,Mk±) 



12e 



5oo(pi - p, M^, Mxo) - Booipi - p, M,, M^) 



- ^(Mi-2M2+p.pi) 

- 3 [M|^o - 2Mlo +P-PI 

- -l^[Ml-2Ml+p-p,) 
+ 6 [p ■ pi - M; 



2 



+ —{Ml-2Ml+p.p,) 



Bu{pi-P,M^o,Mko) 



Bn{pi~p,M^,MKo) 



Bi^{Pi-P,M^±,Mk±) 



5.7 The w-channel contribution 

Finally, the w-channel contribution follows from diagram (e) in figure 01 
1 



5F 



12F2 



{-2Ao{Mko) 



+ 3 



p-p2 + -^{2Ml-2M^ + 5p-p2) 



Bo{p2-p,M^o,Mk±) 



+ -[M^-M'^o-2p-p2 
+ ^{M^^ - Ml - bp . P2) 



+ 2 



+ 3 



3e 



Ml^-p-p2--^{Ml-p-p2) 



Mi„ + — (9Mi - AMI) 



Boip2 - P, M^, Mk±) 

Bo{p2-P,M^±,Mko) 
B^{P2-p,M^o,Mk±) 



+ - [M^ + 2Ml, - AMl± - M^o 



2e 



+ -j={M^^-Ml-bMl) 



B,{p2-p,Mr„MK±) 
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- 4 



5G 



+ 



Mlo-p-p2-^{Ml-Ml) 
8e 



Bi{p2-p,M^±,Mko) 



+ 6(1 + -^) 5oo(p2 - P, M^o, M^±) 



+ 6(^2 + 
+ 10(1 + 



Boq{p2 - P, Mr,, Mk±) 



€ 

^ Boo{p2-P,M^±,Mko) 



- 3 [-M]^o +P-P2 

+ -^(-7Mi + 2M2 + 5p-p2) 

- 3 [-M|-o + 2M^± -P-P2 
+ -^=(Mi + 4M2-5p-p2) 

- 2 [M^± -2M|^o+p-p2 
+ ^{-2Ml + M'^+p.p,) 

^ {-2Ao{Mko) 



Bu{P2-P,M^o,Mk±) 



Bn{p2-p,Mr„MK±) 



Bu{P2-P,M^±,Mko) } , 



12F2 



+ 3 



P-P2 + ^ (2Mi - 2M2 + 5p ■ P2) 
[Ml-Mlo-2p-p2 



Bo{p2-P,M^o,Mk±) 



2e 



+ ^K'-^'-5p-P2) 



+ 2 



+ 3 



3e 



Ml^-p-p2-^{Ml-p-p2) 



Mio + (9M|, - AMI) 
+ ^ [M^ + 2M^o - 4M2± - Ml, 



Bq{p2 -p,Mr„MK±) 

Bo{P2-P:M^±,Mko) 

Bi{p2-p,M^o,Mk±) 



2e 
73 



(M^2 - Mi - 5M2) 



Bi{p2 - p, Mr,, Mk±) 
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- 4 



3e 



M^o-p-p2-^(Mi-M2) 



5i(p2-p,M,±,M;^o) 



6e 

^5oo(p2-p,M,,Mx) 

6 (^1 + -^^ Soo(p2 - P, M^, Mk±) 

2 (^1 + -^^ fioo(P2 - P, M^±,Mko) 



- 3 [-M|o + p ■ P2 

+ -^(-7Mi + 2M2 + 5p-p2) 

- 3 [-M^o + 2Ml± -P-P2 



Bu{P2-P,M^o,Mk±) 



+ —{Ml + AM^-5p-p2) 
- 2 [M^i - 2Mlo +P-P2 
+ ^{-2Ml + Ml + p . p2) 



Bii{p2-p,Mr„MK±) 



Bu{p2 -p,M^±,Mko) 



5.8 Soft virtual photon contribution 

The various topologies of Feynman diagrams containing a virtual photon are 
drawn in figure |3] 

Due to important cancellations between the different contributions from 
these diagrams we will present the result in a compact form, 



6F = j{2Bo{~P2,0,M^)-2Bo{p2+pi,mi,M^) 

- mfCi{-pi,-pi -p^,0,mi,MK) -mfCi{-pi,p2,0,mi,M^) 

- m{C2{-pi,~pi -pu,0,mi,MK) +4^2 ■ piC2{-pi,P2,0,mi, M^) 

- m]{Ml^ - si + 2p ■ pi - 2p ■ p2)Di{-pi,p2, -pi - Pu, 0, mi, M^, Mk) 

- mfiM]^ - si + 2p-pi-2p- P2)D^{-puP2, -pi - Pu, 0, m,, M^, Mk) 

- 2m]pi ■ piDu{-pi,p2, -pi - Pu, 0, mi, M„, Mk) 
+ 2mfpi ■ P2Di2{-Pu P2, -Pi - Pu, 0, mi, M^, Mk) 
+ 2mf{Ml -p-pi+pi-p2-pi-pi) X 

Di3{-Pi,P2, -Pi - Pu,0,mi,M^,MK) 
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+ 2m^pi • P2D23{-pi,P2, -Pi -Pu,0,mi,M^, Mk) 

+ 2m]{M^ -p-pi+pi-p2)D33{-pi,p2,-pi-p^,0,mi,M^,MK)} , 

SG = --{-2Boi-pi,0,mi) 

+ 8p2 ■ PiCo{-pi,p2, m^, mi, M^) + m^Ci{-pi, -Pi - Pu, 0, mi, Mk) 

- {mf - 4:p2 ■ pi)Ci{-pi,p2,0,mi,M^) 

+ m^iC2{-pi,~Pi - Pu,0,mi, Mk) +4^2 ■ piC2{-pi,P2,0,mi, M^) 

- mfiM]^ - si + 2p-pi-2p- p2)Di{-pi,p2, -pi - Py, 0, mi, M^, Mk) 

- m^{MK - si + 2p ■ pi - 2p ■ p2)D2.{-pi,p2,-pi - py,0,mi, M^, Mk) 

- 2m]pi ■ piDii{-pi,p2,~'Pi - Pu,0,mi,M„,MK) 
+ 2mfpi ■ P2Di2{-pi,P2, -Pi - Pu, 0, mi, M^, Mk) 
+ 2m^{M^ - p ■ pi + pi ■ p2 - pi ■ pi) X 

Di3{-pi,P2,-Pi - Pu,0,mi,M^,MK) 

+ 2m^pi ■ p2D23{-pi, P2, -pi - Pu, 0, mi, M^, Mk) 

+ 2mf{M^ - p ■ pi + pi ■ p2)D33{-pi,p2, -pi - py,0,mi, M^, Mk)} ■ 

The infrared divergence in form factors is contained in loop functions with 
m^ in the argument. 



6 Soft photon bremsstrahlung 

Virtual photon corrections to Kj^4 decay rate generate infrared divergencies. 
These cancel, order by order in perturbation theory, with the ones coming 
from real bremsstrahlung corrections. Assume that the emitted photons are 
soft, that is, their energies are smaller than any detector resolution, uj. It 
follows that radiative and non-radiative decays cannot be distinguished ex- 
perimentally and emission of real soft photons should be taken into account. 
Note however that only single soft photon radiation is needed to one-loop 
accuracy. 



6.1 The decay amplitude 

A general feature of photon bremsstrahlung is that, in the soft photon ap- 
proximation, the bremsstrahlung amplitude is proportional to the Born am- 
plitude. Since we deal only with isospin breaking corrections to the F and 
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G form factors, the Born amplitude is taken, all along this section, to be, 

1 



GFV:fi{p,)^,{l-^')v{pi)Q^. 



The contribution of form factors F and G to the Bremsstrahlung amplitude 
can be read off from diagrams in Fig. El 

Let e and q be, respectively, the polarization vector and the momentum of 
the radiated photon. The evaluation of diagrams in Fig. El is straightforward 
and read, to first order in the photon energy, 

^0-7 = eAl- ( ^ - + 0{q) . (103) 

\Prq P2-q J 

Squaring the matrix element (jl03p and summing over polarizations, we 
obtain 



pol. 



m. 



+ 



Ml 



2p2 ■ Pi 



{pi-qY {.P2-qY iP2- q)ipr q)_ 

The preceding expression is singular for vanishing momentum of the soft 
photon. We shall attribute a small but non-vanishing mass to the photon, 
m^, in order to regularize this singularity. 



6.2 The decay rate 

The K£4^ differential decay rate is obtained by squaring the matrix element 
fll03p . summing over spins and polarizations and integrating over the follow- 
ing phase space, 

d% = (27r)^(5(^) iPi+P2+Pi + Pu + q-p) ^ 

d^Pi d^P2 d^pi d^p^ d^q 
(27r)32Ei (27r)32E2 (27r)32E; (27r)32|p^| {2Tif2\q\ ' 

Using the definition of bremsstrahlung integrals as given in the appendix, 
the K(^i^ differential decay rate takes the following form in the soft photon 
approximation, 

spins 
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[I{p2,P2,m^,uj) + I{pi,pi,m^,uj) ~ 21{p2,pi,m^,uj)] , (104) 
where non singular terms have been dropped out. 

6.3 Cancellation of infrared divergencies 

The infrared divergent part of Ke^^ differential decay rate can be extracted 
from (jlU4|) using the definition of bremsstrahlung integrals from the ap- 
pendix, 

"^^^ " 47r2 2Mk ' 

spins 

[l+Pi ■ p2T{-pi,p2,mi, M^)]\nm'^ . (105) 

On the other hand, the infrared divergence coming from virtual photon cor- 
rections to F and G form factors only can be read off from Tab. Q and is 
denoted by dT^^. It is easy then to check that infrared divergencies cancel 
at the level of differential decay rates, 

dr^^ + drf = . 



7 Perspectives 

In this work we studied the decay process, —>■ Tr'^vr^f+z/^, taking into 
account Isospin breaking effects. These come mainly from electroweak inter- 
actions and generate corrections proportional to the fine structure constant, 
a, and to the difference between up and down quark masses, — m^. 

The interest in this decay comes from the fact that the partial wave ex- 
pansion of the corresponding form factors involves vrvr scattering phase shifts. 
The latter can be related in a model-independent way to the vrvr scattering 
lengths which are sensitive to the value of the quark condensate. Thus, a 
precise measurement of form factors should allow accurate determination of 
scattering lengths and, consequently, give precious information about the 
QCD vacuum structure. 

Scattering lengths are strong interaction quantities. On the other hand, 
any Ki4^ decay measurement contains contributions from all possible inter- 
actions, in particular, from electroweak ones. Therefore, it is primordial to 
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have under control Isospin breaking effects in order to disentangle the strong 
interaction contribution from the measured form factors. The present work 
was guided by this motivation and, to this end, analytic expressions for form 
factors were obtained including Isospin breaking effects. These expressions 
are ultraviolet finite, scale independent, but infrared divergent. We showed 
that this divergence cancels out at the differential decay rate level if we take 
into account real soft photon emission. 
Our work should be completed by, 

• a parametrization of form factors in the presence of Isospin breaking, 

• a full treatment of the radiative decay, K^^^. 
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A Loop integrals 

We list here analytical expressions for scalar one-loop integrals and reduction 
of tensor integrals to scalar ones [TH]. Let D be space-time dimension and 'q^^ 
the metric tensor. For large momenta, loop integrals are divergent when D = 
4. A regularization should then be applied to treat this ultraviolet divergence. 
The dimensional regularization consists on calculating integrals for arbitrary 
D. Physical situations are recovered in the limit D —y 4. Infinities show 
up as poles in inverse powers of e = 4 — D > 0. These poles have to be 
absorbed (subtracted) by renormalization constants in order to obtain finite 
(observable) results. This subtraction is not unique and it is a matter of 
taste to chose one of the different schemes. In the MS scheme, one subtracts 
the following. 



A 



- + 1 - 7e + ln(47r) 



(106) 



327r^ 

where 7e is the Euler constant. Another feature of dimensional regularization 
is the scale /z. It has the dimension of a mass and maintains correct units 
while dimensionally regularizing integrals. 



(27r)4 ^ ; (27r 
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Note that the scale // is also absorbed by renormalization constants so that 
obscrvables remain scale independent. 

Some integrals are divergent for vanishing momenta or masses. This kind 
of divergence is called infrared divergence and can be treated in dimensional 
regularization quite as the ultraviolet one {e being negative herein). For 
the situation we are considering, infrared divergence is due to the presence 
of virtual photons. One can then think about a cutoff regularization and 
assigns to the photon a fictitious infinitely small mass, m^, which play the 
role of the cutoff. Obviously, observables do not depend on m^. In fact, the 
infrared divergence due to virtual photons is cancelled by the one coming 
from real soft photon emission leading to cutoff independent quantities. 

In order to evaluate numerically the finite part of loop integrals scalar 
ones are expressed in terms of elementary functions like logarithms. The 
arguments of the latter can always be cast in the following compact form. 



1- Wl- 



4mm' 



(m — m'Y 

a{z,m,m') = \ (108) 



1 + 



bimi' 



z — {m — m') 



l\2 



with z a complex quantity whereas m and m' are real. For complex argu- 
ments, the logarithm is analytic and presents a cut structure given on the 
first Riemann sheet by, 

\n{x + iq) ^ In \x\ + i7re(-a;)sgn(<^) , (109) 

for real x and infinitesimal ^. Herein, is the Heaviside function, 

and sgn is the sign function, 

...J 1 for X > , . 

sgn(<j) = S T r . (Ill) 

[ -1 for X < ^ ' 

For three-point and higher functions, more complicated (but well known) 
functions have to be used. This is the case, say, of dilogarithm, 

U2{z) = - [ dtt-Hn{l- zt) . (112) 
Jo 
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The dilogarithm is analytic in z. That is, it develops an imaginary part for 
values of z fixed by the cut structure of ln(l — z) as seen from the definition. 
From the foregoing one can easily prove the following identity for x real such 
that X > 1, 

Li2(x + z^") — Li2(a; — k) = 2i7r sgn(g) Inx . (113) 

For the convenience of giving compact expressions of loop integrals in the 
various regions of the complex momentum space we will introduce the logical 
function, 

TP/ ,x . f 1 if arrament is true 

liiarqument) = < „ , • r i • (114) 

if argument is false 

The analytic structure of logarithms and dilogarithms is proving relevant for 
the determination of loop functions in physical regions. Let pi and pj be 
external momenta in a given amplitude and define the exchange energy as, 

Pl = {P^-PJ)\ 1 = 1,2,3. (115) 

If mj and nij are internal masses then loop integrals are singular for those 
values given by Landau's equations, 

Pii = {mTrnj)^, (116) 

corresponding respectively to pseudo and normal thresholds of the amplitude. 
Thresholds divide the momentum space into three regions, 

p-j- < (mj - nij)'^ , {mi - mjf < p]- < {uii + nijf , p]^ > {rui + nijf , 

which can be reached from each other by analytic continuation from real to 
complex momentum values, 

P% — ^Pij+^e, (117) 

with e an infinitesimal positive quantity. To see how it works, let us consider 
the analytic continuation of (jlU8|) . 

(Tij = a [p'^j +ie,mi,mj) . (118) 

It is easy to show that the logarithm of the latter function takes the following 
form. 

In [aij) = -If [pI- < [mi - m^Y) x 
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In 



rrii + mj}'^ - pj- - ^/ {rrii - mjY - pf^ 



+ 2i If [{rrii - 'rrijY < Pij < ("^i + ^jY) x 



arctan 



Pij - {rrii -TTijY 



\J i^i + mjy - p% 



- \i{pl> {rrii + mjf) 



In 



\Jp% - {mi-mjY - ^pj. - {nii + nij) 



.^2 



;ii9) 



Finally, the following notations are usefuU for the reduction of vector and 
tensor integrals to scalar ones [7^] , 



= f -ml + ie, Ni = {pi + if - rrf + ie , 



(120) 



2p,-l = N,-No- f,, fi = pI -m1 + ml. 



;i21) 



A.l One-point functions 

The one-point functions or tadpole integrals are defined by, 



4-D 



{27t)d 

In dimensional regularization the scalar integral reads 

1 



An = rrin 



-2A 



The vector and tensor one-point functions are found to be. 



rrin 



Tfl 

2Ao(mo) + — ^ 



(122) 



(123) 



(124) 
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A. 2 Two-point functions 

The two-point functions are defined by, 

{5o,5^5'^^}(pi,mo,mi) = 

The scalar integral reads [HUj , 

_ 1 Ao(mo) 1 Ao(mi) 

-^0 ~ o 2 o 2 



+ ~ V ("^0 + mi)'^ - Pi V ("^0 - mi)^ -pi y< 
A/(mo + mi)^ - + A/(mo - mi)2 - pj 

ill . = . = X 

(mo + mi)2 - p2 _ ^(^^^ _ _ p2 

If (pi < (mo - mi)^) 
2 



2 ^ (mo + mi)2 - pi \JpI- (mo - mi)^ x 



Pi 

arctan ^^^^^^2 x 
y/(mo + mi)2 - pi 

If ((mo - mi)^ <p\< (m-o + mi)^) 

^ \JpI- {niQ + mi)2 y^p? - (mo - mi)^ x 

a/pI - (mo - mi)2 + y/pl - (mo + mi)^ . 

m — — — tTT 

^y pI - {mo - rriiy - ^y pI - {mo + 

If (p^ > {mQ + mif) 



The vector two-point function is written as, 

= p^B,, 

with the coefficient, 

2plBi = Ao{mo) - Ao{mi) - f^Bo 
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The tensor two-point function possesses the following decomposition, 

5/^- = r7^-i?oo + KKi?ii, (129) 

with the coefficients, 

18^00 = 3Ao(mi) 

167r2 

^2 



+ QmlBo + 3/i5i + -i^ {Ami + 2^ - /i) , (130) 



ISpfBn = 6Ao(mi) 

- QmlBo - l2hB, - -1^ [Ami + 2m? - h) . (131) 

Two-point functions are ultraviolet divergent. It is convenient to define ultra- 
violet finite parts by subtracting the corresponding poles. The finite parts are 
noted with the superscript "r" like renormalized and are defined as follows, 

5o = S5-2A, Si = 5[-A, (132) 

Eoo = Bl^ + ^ (P? - 3m2 - 3m2) A , B^ = Bl^-\\. (133) 

A. 3 Three-point functions 

The three-point functions are defined by, 

{Co,C^,C^''} (pi,P2,"^0,"^l,"^2) = 

(0j(^oiViiV2)-Ml,^^^"^'^}- (134) 

The general expression for the scalar integral is complicated especially 
for separating real and imaginary parts. In order to make this separation 
more transparent we will give a one-dimensional integral representation of 
the scalar integral 



16vr^ Jo 

I T-/ 2 2 2 2 2 2\ 

+T{x; P2, Pi,Pi2, rriQ, m^, m^) 
+J^{x;pl2,pl,pl,ml,ml,ml)] . (135) 
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Herein, 



with, 



4a 



n~t 2 2 2 2 2 2\ ■ 



ln(2a + & + d) - ln(6 + d) 
+ If(A < 0) 1^ Inc- ^ ln(a + 6 + c) 



+ 



— arctan 



arctan 



+ If(A>0) 



7^ 
2a + 6 



+ 7rlf(6 < < 2a + 6) 
ln|2a + 6- v^A] 



In |6 - v^l - m\i{h < < 2a + 6)] 



2VA 



ln|2a + 6 + VAI 



ln|6 + VA| +i7rlf(6 < - VA < 2a + 6)j } , 



a 
6 



/'2i 2 2\i 2, 2 2 

(m^ + m2 - V\2)^ + ^2 + mo - . 



2 2 I / 2 I 2 2\ I 2 

c = m^x + [rrii + uIq — p^)x + , 
d = 2a(l + x)-6, 
A = 6^ - 4ac + ie . 



(136) 



(137) 
(138) 
(139) 
(140) 
(141) 



Obviously, this representation is appropriate for direct numerical integration. 
The vector three-point function is written as. 



The coefficients are given in a compact form. 



a = E(^2)rnC„, i = 1,2, 



(142) 



(143) 



n=l 
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where, 

2pl 2pi ■ p2 



2pi ■ p2 2p^ 



(144) 



and, 



Ci = Bo{p2,mo,m2) - Bo{p2-pi,mi,m2) - fiCo, (145) 
C2 = Bo{pi,mo,mi) - Bo{p2-pi,mi,m2) - f2Co. (146) 

The tensor three-point function possesses the following decomposition, 

2 

C"^^ = r^Coo + E P^PlCj, . (147) 
j,k=i 

The coefficients are found to be, 

4Coo = BQ{p2-pi,mi,m2) + 2mlCo + fiCi + f2C2 + Y^, (148) 
2 

C,,- = 5^ (^2)7,' (Q„ - 25,„Coo) , ^,J = 1, 2, (149) 

n=l 

where, 

Cii = Bo{p2-pi,mi,m2) + Bi{p2-pi,mi,m2) - fiCi, (150) 

C12 = Bi{pi,mo,mi) ~ Bi{pi-p2,m2,mi) - f2Ci, (151) 

C21 = Bi{p2,mo,m2) - Bi{p2- pi,mi,m2) - fiC2, (152) 

C22 = Bo{pi-p2,m2,mi) + Bi{pi-p2,m2,mi)- f2C2. (153) 

Only the tensor coefficient Cqo is ultraviolet divergent. The finite part is 
defined by the following subtraction, 

Coo = Co^o-^A. (154) 

A. 4 Four-point functions 

The four-point functions are defined by, 

{Do,D'',D^''} {pi,p2,P3,mo,mi,m2,m3) = 



46 



-i^'-"" I ^{NoN.N^N.r'iU^J^n . (155) 



The expression for the scalar integral is known but difficult to handle 
Since we are interested only in the case, mo = rrij, pi = ml, p\ = m\, 
the general expression will not be reproduced here. The vector four-point 
function is written as, 

= p>tDi+ D2 +p'^D3. (156) 

The coefficients are given in a compact form, 

3 

A = $^(23)r„'^^n, 2 = 1,23, (157) 



n=l 

where, 

/ 2pl 2pi ■ p2 2pi ■ p3 
Z3 = 2p^-p2 2pl 2p2-p3 I , (158) 
\ 2pi ■ p3 2p2 ■ P3 2pl 

and, 

T^i = Co{p2,P3,mQ,m2,m3) 

- Co{p2-pi,P3-Pi,mi,m2,m3) - fiDo, (159) 
T>2 = Co{pi,p3,mo,mi,m3) 

- Co{p2 - Pi,P3 - Pi,mi,m2,m3) - f2Do , (160) 
1^3 = Co{pi,p2,mo,mi,m2) 

- Co{p2-pi,P3-pi,mi,m2,m3) - f3Do. (161) 

The tensor four-point function possesses the following decomposition, 

3 

D'^'' = r^'^'^Doo + J2 P'PkDjk . (162) 
j,k=i 

The coefficients are found to be, 

3 

2Doo = Co{p2-pi,P3-Pi,mi,m2,m3) + 2mlDo + '^fnDn, (163) 

n=l 
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3 

Dij = J2 (^3)7i (An - 25,,Doo) , ^, J = 1, 2 3 , (164) 

n=l 

where, 



2^11 




Co 


\jr ^ Jr ^7 Jro ^ i.^ • • ^ l 


mo, 1713) 


+ Ci(p2 


— V^,V3 


— 'P^,m^ 


mo, 777^) 




1 


Co 


Tin — T)i T)i — T)i 777,1 


777,0 777, q1 


- f^D^ 
J 1^1 -1 






(165) 






Ci 


Pi , Ps, mn, 777i , 777s) 
















Co 




777,0 777,1 1 








(166) 


2^13 




Ci 


{pi,P2,mo,mi,m2) 
















Ci 


[Pl -P3,P2-P3,m3 


nil, m2) 


J 3-^1 ) 






(167) 


i^21 




'-'1 


{P2:P3: mo, 1712, 1713) 
















Ci 


[P2 -pi,P3-Pi,mi 


1112,1713) 


-flD2, 






(168) 


2^22 




Co 


{P3 -P2,Pl -P2,m2 


1713,1711) 


+ Ci(p3 


-P2,Pl 


-P2,m2 


7773,7771) 




+ 


C2 


{P3 -P2,Pl -P2,m2 


1713,1711) 


-/2i^2, 






(169) 


^^23 




C2 


{pi,P2,mo,mi,m2) 
















C2 


[Pl -P3,P2-P3,m3 


mi, 777,2) 


-f3D2, 






(170) 






C2 


{P2,P3,mo,m2,m3) 
















C2 


[P2 -pi,P3-Pi,mi 


1712,1713) 


-hD3, 






(171) 


2^32 




C2 


[pi,P3,mo,mi,m3) 
















Cl 


{P3 - P2,Pl - P2, ^2, 1713, mi) 


- f2D3 , 






(172) 






Co 


[Pl -P3,P2 ~P3,m3 


1711,1712) 


+ Ciipi 


-P3,P2 


- P3: m3 


7771,7772) 




+ 


C2 


[Pl -P3,P2-p3,m3 


1711,1112) 


- f3D3 . 






(173) 



A. 5 Infrared divergent integrals 

From all integrals cited before only the scalar three- and four-point integrals 
are infrared divergent for the particular case, 

777o = 777^ , Pl = 777i , p^ — m^ . (1^4) 

Furthermore, these integrals are related by, 

Co(pi,P2,m^, 777i,7772) = lim [-7773i:>o(pi,P2,P3,"7T,, 7771, 7772, 7773)] ■ (175) 
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In this subsection, the integrals in question will be denoted Cq and D^, 
respectively, and read jHlI , 



Dn 



1 



1 



0"12 



IGtt^ mim2 1 — C7^2 
In (cri2 



X 



1 



21n(l-ay--ln(ai2) 



In 



f ^1 \ 
\m1m2 J 



+Li2 ( 1 - 0-12 



nil 

777,2 
1 



Li2 

0"12 



1 - (T12 



772,2 



:i76) 



167r^ 77^l77^2 p| — ml 1 — cr^2 



2 In ((712) 
2 



In (1 - a^^) - In 



777 



3 - - «e 



+ - 



TT 



Li2 ((Ti^2) + (^13) + In^ (^23) - [Lia (CTi2af3a^3) 



k i 



+ (In (ai2) + In {a',,) + In {ai,)) In (l - (712^^3(7^3)] } . (177) 

The In 777^ terms in the preceding formulae generate infrared divergence in 
the virtual photon correction to the amplitude. This divergence is cancelled 
by the one coming from the associated soft real photon emission originating 
from the bremsstrahlung integrals. 



2pi ■ P2 



(27r)32|g| (pi-g)(p2-g) 



(178) 



The expression of bremsstrahlung integrals in terms of logarithms and dilog- 
arithms reads, 



api ■ p2 



Atx"^ api ■ P2 — P2 



2api ■P2~pI 2c^ 



p?- 


\P1 


1 


Pl + 


\P1 


~ 4 



.Pi 



\P2 



P2 + \P2 



Lio 



-Li2 (l-^(p^+|P2l)) 
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Li2 (l-^(p?-|Pil))-Li, 



(179) 



where, 



P = T ~f' ■ (180) 



and a is defined through, 



aYi - 2api -^2+^2, n > • (181) 

Pi 

For the particular case, pi = p2 the integral simplifies to, 

/(pi,pi,m„a;) = ^ An ^ + ^ In 4^^) . (182) 

In order to make explicit the cancellation of infrared divergence in the am- 
plitude we will use a compact unified notation for the Inm^ terms in the 
integrals Cq, Dq and /. To this end, define the r function to be, 

/■I dx 
t{pi,P2, nil, 1712) = / -^-^ r^— 2 2^ — \ 2- (1^3) 

The r function can be expressed in terms of the a function as follows, 

r{pi,P2, 1711,1712) = — ^ ln((Ti2) , (184) 



11711712 1 — a 



where the logarithm can be read from (jll9|) . In terms of the r function, 
the infrared divergence in the three-point function, four-point function, and 
bremsstrahlung integral reads. 



Co — > T {pi,p2, nil, 1712) In in^, (185) 



-Do — ^ ^2 J ^2 ^(Pi^P2,mi,m2) Inm^, (186) 

/ — > -^pi-p2T{pi,p2,nii,ni2)lnm^^. (187) 



327r^ pI — m 



3 
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Figure 1: Angles and unit vectors used in the kinematical description of 
decays. Etttt and E;,^; are the planes defined in the kaon rest frame by 
the pion pair and the lepton pair, respectively. 9^^ (Oi), the angle formed 
by Pi (Pi), in the dipion (dilepton) rest frame, and the line of flight of the 
dipion (dilepton) as defined in the kaon rest frame. 0, the angle between the 
normals to E^^ and Y^iui- v is a unit vector along the direction of fiight of 
the dipion in the kaon rest frame, c (d) is a unit vector along the projection 
of Pi (pi) perpendicular to v. 



diagram 5G 
3. (a) ^^"^5 



Table 1 : Infrared divergent part of the corrected / and g form factors due to 
virtual photon corrections. The contribution from diagram 3. (a) comes 
from wave function renormalization of external charged particles, tt", and 

e+. 
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(d) 



Figure 2: Feynman diagrams representing decay amplitudes at tree level. 
Wavy lines stand for photons. Only diagrams (a) and (b) contribute to the 
decay amplitudes A^^ and . 
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1+ 1+ 



(d) (e) (f) 

Figure 3: Feynman diagrams representing the Kf^ decay amplitude of the 
neutral kaon at one-loop. Is shown only the non photonic topology. Diagram 
(a) represents Born and counter-terms contributions. Diagram (b) accounts 
for 7r° — 77 mixing. Tadpole contribution is given by diagram (c). Diagrams 
(d), (e) and (f) stand for contributions from the s-, t- and li-channels, re- 
spectively. 
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TT '^l TT 




(d) (e) (f) 

Figure 4: Virtual photons in Feynman diagrams for the Ku decay in the 
mixed channel. Diagrams (a) and (b) represent tadpoles. Diagram (c) gener- 
ates three-point functions of type leg-leg photon exchange. Diagrams (d) and 
(e) generate three-point functions of type vertex-leg. Four-point functions of 
type leg-leg are generated from diagram (f). 



61 




Figure 5: Feynman diagrams representing the contribution of F and G form 
factors to the bremsstrahlung amphtude. 
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